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Abstract 

We show that by adding a vector-like 5 + 5 pair of matter fields to the spectrum 
of the minimal renormalizable SUSY SU(5) theory the wrong relations for fermion 
masses can be corrected, while being predictive and consistent with proton lifetime 
limits. Threshold correction from the vector-like fields improves unification of gauge 
couplings compared to the minimal model. It is found that for supersymmetric 
spectra lighter than 3 TeV, which would be testable at the LHC, at least some of 
the nucleon decay modes should have partial lifetimes shorter than about 2 x 10 34 
yrs., which is within reach of ongoing and proposed experiments. 
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1 Introduction 



While elegant and simple, the minimal renormalizable supersymmetric S77(5) model [1, 
2, 3] suffers from two main drawbacks. The first is the wrong predictions it makes for the 
light fermion masses. This theory predicts the asymptotic relations m d = m°, = m[J 
and m° = m° connecting the charge —1/3 quark masses and charged lepton masses, valid 
at the grand unification scale of 2 x 10 16 GeV. Such relations would enable one to calculate 
the down-type quark masses in terms of the charged lepton masses by evolving the mass 
parameters via the renormalization group equations (RGE). The relation = mP T is 
generally considered a successful prediction of minimal SUSY SU(5), since the 6-quark 
mass computed in terms of r-lepton mass is typically within about 20% of its experimental 
value. The relations involving the lighter families, however, lead to wrong predictions. For 
example, the RGE-invariant relation m,d/m s = me/m^, which follows from the asymptotic 
relations of the minimal model, differs from experimental values by about a factor of 10 
(md/m s ~ 1/20 while m e /m M ~ 1/200 at low energy scale [4]). 

The second drawback of the minimal SUSY SU (5) model is its prediction for proton 
lifetime for the mode p — > VK + which arises via the exchange of colored Higgsinos. The 
lifetime is generically too fast compared to the present experimental limits. This predic- 
tion follows mainly from the requirement of gauge coupling unification. The spectrum of 
the minimal supersymmetric standard model (MSSM) at low energies does not lead to a 
precise unification of the three gauge couplings when the full two-loop RGE are used, and 
therefore requires some threshold correction from the GUT scale. The only possibility in 
the minimal renormalizable SU (5) set-up is to make the color triplets from the 5# + 5# 
Higgs fields (which transforms as (3,1,-1/3) + h.c. under SU(3) C x SU(2) L x U(1) Y 
gauge group) somewhat lighter compared to the vector supermultiplets (the X and Y 
gauge bosons of 577(5)). Since the same color triplets mediate d = 5 proton decay [5, 6], 
making it lighter than the GUT scale results in a considerably shorter proton lifetime 
[7, 8, 9, 10, 11, 12], typically in conflict with experimental limits. Notice that this out- 
come is due to the minimal particle content: the same color triplet that corrects the RGE 
running of the gauge couplings is coupled to the Standard Model (SM) fermions with 
fixed Yukawa couplings. (The color triplet Yukawa couplings are unified with the Yukawa 
couplings of the SU(2) L doublets also contained in 5# + 5# that generate quark and 
lepton masses and mixings.) There is no other choice in the minimal model for correcting 
the RGE running of the gauge couplings. 

There are various well known ways out of these two problems. The most commonly 
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used solution is the inclusion of higher dimensional operators. Due to the vicinity of 
Mqut to M Planck such operators may not be negligible numerically, especially for the 
lighter fermion masses [13]. For example, they can easily improve the calculated masses 
of the first two generations. Their influence for proton decay is even bigger. They make 
the Yukawa couplings to the color triplet Higgs different from those to the weak doublet 
Higgs, so that there is some freedom which can be used to somewhat suppress the d — 5 
proton decay amplitudes. Alternatively, these higher dimensional operators can allow for 
a lighter color octet and weak triplet (remnants of SU (5) symmetry breaking via a 2A H ) 
which can increase both the GUT scale and the color triplet masses [14, 15, 16], alleviating 
the d — 5 proton decay problem significantly. 

The problem with this natural solution is that it automatically introduces a large 
number of new parameters into the game, thus precluding any quantitative prediction. 
So, although the model can be made consistent and realistic, it is difficult to test it. 
There is also some questions about the strengths of these higher dimensional operators 
being of the right magnitude if they are induced by quantum gravity effects. In this 
paper we take a different approach. We assume that our supersymmetric SU(5) GUT 
is renormalizable. After all, we really do not know how gravity influences our particle 
physics world, and a conservative approach would be to not rely heavily on gravity- 
induced corrections. This approach of using only renormalizable couplings has brought 
great success in the electroweak sector of the Standard Model. The renormalizability of 
the theory would greatly reduce possible couplings in the theory resulting in enhanced 
predictivity. With this in mind we shall add to the minimal supersymmetric SU (5) as little 
as possible: a vector-like 5 + 5 matter field. This will allow unequal mixings of the down 
quarks and charged leptons with these fields, thus correcting the wrong mass relations. 
Simultaneously this set-up would provide a new set of color triplet /weak doublet fields, 
which allows for a precise unification of gauge couplings by choosing the color triplet 
somewhat lighter than the weak doublet. Note that such a choice does not run afoul with 
d = 5 proton decay rates, unlike the minimal SUSY SU (5) model, since the 5 + 5 fields do 
not acquire vacuum expectation values (VEVs). As in minimal SUSY SU(5) we assume 
.R-parity conservation, and we take the vector-like 5 + 5 pair to be fermion-like. Had we 
chosen Higgs-like multiplets such as 45 + 45, the wrong fermion mass relations could have 
been corrected [17], however in this case quantitative predictions for proton decay would 
be difficult to make owing to the large number of parameters that would be introduced. 
Another possible solution to the wrong mass problem of the minimal SUSY SU(5) model 
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is through supersymmetric threshold corrections arising from soft SUSY breaking terms 
with a particular form, see for example Ref. [18, 19]. Here we shall assume that the SUSY 
spectrum is such that such threshold corrections remain small. Yet another possibility 
is to utilize large Yukawa couplings involving vector-like multiplets. This can raise the 
unification scale when two-loop RGE effects are included, which would allow for a better 
prediction for a 3 (M z ) [20, 21]. 

We now turn to the discussion of fermion masses in presence of a 5 + 5 matter fields 
and show how the mixing of these fields with the MSSM fermions corrects the wrong mass 
relations. We then derive the baryon number violating effective d — 5 superpotential and 
study its implications for nucleon lifetime. The small number of new parameters that 
are introduced with the addition of a 5 + 5 vector-like fermions allows the model to be 
consistent with current proton lifetime limits, but at the same time we find that at least 
some modes should have partial lifetime less than about 2 x 10 34 yrs. In our analysis we 
assume that the GUT scale stays well below the Planck scale (by a factor of 20 to 50) 
so that quantum gravity effects can be ignored, and the approximate unification of the 
gauge couplings that occurs in the MSSM is not a complete accident. For supersymmetric 
spectrum, we assume that all super-particles have masses less than about 3 TeV, which 
would make them detectable at the LHC, while at the same time providing a solution to 
the gauge hierarchy problem. 

2 Fermion masses with vector— like 5 + 5 matter fields 

Before discussing the modifications of the fermion mass relations with the inclusion of a 
5 + 5 matter fields in SUSY SU(5), let us briefly summarize the situation in the minimal 
renormalizable SUSY SU (5) model. 

2.1 Fermion Masses in minimal SUSY SU(5) 

The matter fields of the model consist of three generations in representations 10j + \, 
i = 1,2, 3. The Higgs sector consists of an adjoint 2A H used for breaking £77(5) symmetry 
down to the SM symmetry, and a pair of 5# + 5# fields for electroweak symmetry break- 
ing. The renormalizable superpotential of the adjoint field relevant for SU (5) symmetry 
breaking is 

W 24 = ^Tr (24^) + ±Tr (24%) . (2.1) 
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The scalar potential induced by this superpotential has a ground state with a non-zero 
vacuum expectation value, 



(24 H )=t;diag(2,2,2,-3,3) (2.2) 

which spontaneously breaks S77(5) -^■SU(3) c xSU(2) i xU(l)y. The VEV v is determined 
to be 

(2.3) 

The simplicity of Eq. (2.1) fixes the masses of the color octet (the (8, 1,0) fragment of 
24h which is a physical Higgs particle) M 8 and the weak triplet (the (1,3,0) fragment of 
2A H ) M 3 to be 

M 3 = M 8 = 5m . (2.4) 
The same VEV sets the super-heavy 377(5) gauge boson masses to be 

M x = My = 5^2^ . (2.5) 
A 

The two MSSM Higgs doublets H u and live in the pair of Higgs fundamentals 
5/f + 5h and have Yukawa couplings with the matter fields given by 

W Y = lO^lO^ + biYi^jbH . (2.6) 

The equality of the down-type quark masses and charged lepton masses follows from this 
superpotential: 

M D = (5 H )Y? = Ml . (2.7) 

The color triplets from 5 H + 5# have the same Yukawa couplings as the Higgs doublets 
and would mediate rapid proton decay via d = 5 baryon number violating operators. For 
this reason they must be ultra-heavy, preferably with a mass above the GUT scale. In 
the superpotential terms 

W 5 = 5 H (m H + Vh24 h ) $h (2.8) 
this can be arranged by a fine-tuning: 



The color triplet mass is thus 



7TI 

m H = 3r] H j . (2.9) 



7TI 

M T = 5ri H — (2.10) 
A 
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which shows that m T cannot be arbitrarily large if we demand (as we do) perturbativity 
of the couplings: 

Due to the relation in Eq. (2.4), the requirement of gauge coupling unification would 
imply that the color triplet mass is actually much lower, around or even smaller than 10 15 
GeV [ll]. 4 Such a light color triplet would mediate too fast a proton decay, which is a 
problem with the minimal model. 

2.2 Mixing of chiral families with 5 + 5 fields 

To the minimal SUSY SU(5) described in the previous subsection we now add a vector- 
like pair of matter fields 5 denoted as 5 4 + 5 4 . With their i?-parity assumed to be identical 
to that of the chiral families 10j + 5j (or equivalently odd matter parity), the most general 
renormalizable addition to the superpotential of minimal SU (5) is 

W 4 = 5 a (pt a + Va 2A H ) 5 4 , a = 1, . . . , 4 . (2.12) 

Notice that, without loss of generality, by an appropriate choice of the basis, the terms 
5 4 10j5# can be rotated away. Thus, the whole Yukawa superpotential reads as 

W Y = 10^10,5^ + o,V-; ; l ()/,„ + 5 a {n a + r] a 2A H ) 5 4 . (2.13) 

One can work in a basis where the 3x3 coupling matrix Y^ 3 is diagonal: 

Y5 3 = yAj ■ 

Plugging the VEVs (5 H ) = v u , (5 H > = v d ,(24 H > = v diag (2,2,2,-3,-3) into Eq. 
(2.13) and keeping color triplet states T, f (from 5h,5h), the relevant terms involving 
the MSSM fields and the additional vector-like states will be 

W Y = L T Mf xi E c + D cT Mf^D + « T M°« C + l T Y 5 qf + —u T M°dT 

+ d cT Y 5 u c T + — e cT M?jU c T , (2.14) 

V u 

where 

L T = (h,h, h, h) , E cT = {e\, e c 2) e%, U) , 



4 An exception would be to choose very special MSSM soft parameters [22]. This may however require 
very particular and exotic hidden and messenger sectors of SUSY breaking. 

5 The use of heavy vector-like matter to correct the bad mass relations in GUTs is long known. For 
an incomplete list see for example [23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34]. 
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D cT = (dl, d c 2 , d c „ d%) , D T = (d ± , d 2 , d 3 , d\) , 



(2.15) 



M 4x4 = / VAjVd M\ \ M 4 X4 = ( ViSijVd M? \ (2.16) 



o |Afii ; ' d v i-4i 

M\ = n a - 3 Vi v , Mf = fjn + 2 Vi v , M° = F 10 t; u . (2.17) 

Let us now focus on the light (MSSM) charged lepton and down-type quark masses 
arising from Eq. (2.16). These are obtained by removing the heavy vector-like state from 
the spectrum. The mass matrices of Eq. (2.16) can be block-diagonalized so as to bring 
the mass terms in the superpontential to the form 

W mass = e T M E e c + d T M D d c + u T M v u c + M D DD + M C CC . (2. 18) 

The reduced mass matrices M E and M D , derived in Appendix A.l, can be made real and 
have forms 



d x c\ \ / d lC f -d lS d s d 2 -d^sfs 

M E = I -d^sl d 2 c e 2 , M D = \ d 2 c d 2 -d 2 s 2 l s d 
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-d\c%s\s% — d 2 s 2 sl d^c\ I \ d 3 c d 
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(2.19) 



with 



\ym\ 



, ^ = cos^, s^sin^, tf^tanff 



j.e,d 1-^1 ,e,d I -^2 I „e,d ,e,d Jj^3_[ e,d e,d rn nr.\ 

1 " \M l /\ ' 2 " \M 1 / 1 1 ' 3 " \M l /\ 2 ■ { ] 



Note that since M\ ^ M d , the wrong GUT scale asymptotic relation M E (M G ) = M%(M G ), 
which is problematic for the minimal renormalizable SU(5) model, is avoided here. In 
Eq. (2.16) Mu = Mjj = Y w v u , since the up-type quarks do not mix with any of the 
vector-like field. 

From Eq. (2.19), it follows that realizing the mass hierarchy between different families 
is possible only when the diagonal factors di are hierarchical, d\ <C d 2 <C d 3 , in which case 
we can write down very simple formulas for the masses: 

ml* ~<U cos % ,d ■ (2.21) 

Thus, it is possible to fit all quark and lepton masses consistently to the observed values. 
The mixing angles are related by the ratios: 



The 3x3 light fermion mass matrices are diagonalized via bi-unitary transformations 
M E = VlM*V E , M D = tf D M°V D , Mu = VX aq V:, (2.23) 



diag E i 1V1 D D diag 

by going from the flavor to the mass eigenstate basis: 6 

d -> UlPd , e -> E/f e , u ->■ VfP 1/2 u , 1/ -> E/f 1/ 

cf ->■ V^P*cf , e c ->• V ] E e c , m c ->• VjVP*u c . (2.24) 

The diagonal phase matrices P and P are introduced (see Appendix A.l for details) so 
that the CKM matrix can be written as 

Vckm = y/P~*V:ulP (2.25) 
in a standard parametrization with a single phase: 



C12C13 S12C13 s 13 

VcifM = I — s 12 c 23 ~~ c 12 s 23^13 c 12 c 23 ~~ s 12 s 23Sl3 S23C 13 

S12S23 - C12C23S13 - C12S23 - S12C23S13 C23C13 



(2.26) 



The entries of Eq. (2.26) can be parameterized by four Wolfenstein parameters A, A, p 
and fj as follows: 



sn = A , ci 2 = Vl- A 2 , s 23 = ^A 2 , c 23 = Vl - A 2 A 4 

With the central values of these parameters taken from PDG [37] 

A = 0.2253, A = 0.808, p = 0.132 , ?? = 0.341 (2.28) 

we can calculate the CKM elements at M z scale. The corresponding CKM elements at 
the GUT scale are obtained from Vckm(M z ) by dividing the 13, 23, 31 and 32 elements 
by a common RGE factor (~ 1.055 for tan (3 = 7), while keeping the remaining elements 
intact. 



6 Neutrino masses are ignored for simplicity, since they are irrelevant for our studies. They can of 
course be included via the seesaw mechanism with right-handed singlet neutrinos fields introduced. This 
would have very little effects on our discussions. Another possibility would be to include bilinear R-parity 
violating couplings, see for example [35]. 
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As far as the charged fermion masses are concerned, their Yukawa couplings at the 
GUT scale, taken to be M G m 2 ■ 10 16 GeV, for tan/3 = 7, are taken to be 

Kiag/vu = diag (5.49 ■ 1CT 6 , 0.00323, l) X t (A G ), X t (M G ) ~ 0.44 , 

M? iag /v d = diag (0.000886 , 0.01646 , l)X b (M G ), X b (M G ) ~ 0.038 , (2.29) 

^Ia S /^ = diag (0.0002777, 0.05862, 1) X T (M G ), X T (M G ) ~ 0.047 . 

These values correspond to central values of these masses at low energy scale, see for eg., 
Ref. [36]. These numerical values will be used below for the study of proton decay. We 
emphasize that realistic fermion masses are obtained in this model, unlike the minimal 
renormalizable 577(5) model. 

3 The value of a 3 (M z ) 

Since in the model under study we have additional states D,D,C,C beyond those of 
minimal SUSY S77(5), if their masses lie below the GUT scale (M G ), the unification of 
three gauge couplings will be modified. The masses of these extra states are given by 



M D = \J\M[\ 2 + \M l 2 \ 2 + |M|| 2 + \M\\ 2 , 



M c = sJ\Mf\ 2 + |M|| 2 + |M|| 2 + \Mf\ 2 . (3.1) 



since m M l a , M d a there are S77(5) symmetry breaking effects (see Eq. (2.17)), in general 
these two masses differ: Mp 7^ M G . We will exploit this fact for improving the value of 
a^{Mz) predicted by the demand that the three gauge couplings unify. Assuming that 
M D ~ M G and M c < M G , we will have: 

1 - 1 9 >n^, (3-2) 



a 3 {M z ) al(M z ) 14tt M g 

where a°(M z ) denotes the value of the strong coupling constant one would have obtained 
in minimal SUSY SU (5) GUT. The second term on the right-hand side of Eq. (3.2) is due 
to the one-loop contribution of the extra color triplet pair from the vector-like fermions 
with mass M G < M G . With the choice of super-particle spectrum inspires by supergravity 
(see below Eq. (4.20) and Table 1 for the spectral values we use), and with all the GUT- 
scale states (besides (7,(7) having masses ~ M G one would obtain a\{M z ) ^ 0.127. To 
bring this somewhat large value down we take ^ ~ 0.061. Using this in Eq. (3.2), 
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we obtain a 3 (M z ) ^ 0.1184 - the central value of the experimentally determined strong 
coupling constant. 

Note that from Eq. (3.2) the ratio ||g is determined. The value of M G should be found 
from the meeting point of three gauge couplings. Because of the fact that the dependance 
of M G on a.i(M z ) is exponential, we are able to determine M G) and therefore also M T , 
only to an accuracy of about 22%. This will cause an uncertainty of about 45% in the 
d — 5 proton decay lifetime estimate. Further uncertainty is caused by the uncertainty 
in the ratio r = M 8 /M X - The natural value of r is of order one, but r C 1 cannot be 
excluded. Choosing r < 1 would result in larger values of the unification scale, which 
we shall demand to lie at least a factor 20 — 50 below the Planck scale, so that quantum 
gravitational corrections to the gauge coupling evolution remain small. 

4 Effective baryon number violating operators and 
nucleon decay 

In studying nucleon decay, we will need to derive the relevant d — 5 baryon number 
violating effective operators. These operators are obtained by integrating out the extra 
vector-like matter superfields, as well as the states T, T from the couplings given in Eq. 
(2.14). Details of this procedure are given in Appendix A. 2. Here we present the relevant 
effective superpotential couplings: 

W eff = W mass + Wt 5 + Wt 5 , (4.1) 
where W mass is given in Eq. (2.18), 



,abc 

Wt 5 = ~T7 (u T a Mud b )(is T M E P'd c - e T M E P'u c ), (4.2) 

M T v u v d 



and 

f abc 

W d R = 5 = T7 (ufMuP'*e c )(dfMEK) . (4.3) 

M T v u v d 

Here a,b,c are color indices. P' is a phase matrix P' = diag(e l<Sl , e^ 2 , 1). M D and M c 
are the masses of the extra vector-like weak doublets (D, D) and color triplets (C, C) 
respectively. Note that all these coupling are written in the flavor basis of MSSM quarks 
and leptons. 7 The couplings given in (4.1)-(4.3) will be needed for the discussion of 
nucleon decay. Now we turn to the estimate of d — 5 proton decay rates. 



7 These states differ from those of initial superpotential (2.14) due to various rotations (discussed in 
the Appendix). However, in Eqs. (4.1)-(4.3) we use the same notation (without primes) for simplicity. 
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4.1 Effective d = 5 operators in the mass eigenstate basis 

With the basis change given in Eq. (2.24) and using Eqs. (2.23), (2.25), the baryon 
number violating operators of Eqs. (4.2), (4.3) will have the following form in the mass 
eigenstate basis: 

~abc , v 

< =5 = {u T a PM u dtag V CKM d b ) [v T M E diag Vd c - e T Mi ag VVi KM u c ) (4.4) 

,abc , . 

= j^- d {uf M u diag V C KMV^e c ) (d?M° ag Vl KM P*ul) . (4.5) 

The matrices V and P are given in Eqs. (A.10)-(A.12). 

The d — 6 four fermion operator obtained from Wl =5 by wino dressing and involving 
the neutrino has the form 

f abc 

Ott = —r 5aip {u 8 a dt) (rfX) , (4.6) 

where 

CLr, = 92 E - c ^p) I h=Mq (V C KM)p a (V£ KM ) Sa I(u?, cT, W)A s (d\ J, d°) 

+9 2 2 &*Pp ~ c^p)\»=Mc ( V ckm)^ I(u p , e p , W)A s (d a , u\ v?) , 
, c Palp = -i- {M u diagPVcKM )^ (V T M° ag ) ip , 



V u V d 

Qfctf, = -i- (Ml/M fo (^M^ T Mf a J . (4.7) 



Here / is the loop integral defined as 

I(i 3 k) = _i ( m " l„ !^ _ m ? l n (4.8) 

' ' 167r 2 m\ - m 2 - \mf - m\ m\ m? - m\ m 2 k J ' 

while A# accounts for short distance renormalization factor of the corresponding LLLL 
d=5 operator. Here we present some of these RG factors, which will be needed later on 
for numerical calculations: 

A s (d\uP,d a )^^ 3 = A s (d a ,u s ,uP) aA ^ 3 ^ 6.88 , 

A s {d? -y,b) y , p ^ = A 5 (cF,t,cT)w3 = A s (d a ,u s ,t) a ^ ^ 6.54 , 

A s (d\t,b)^ 3 ~ 6.2. (4.9) 
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These expressions are valid for low to moderate values of tan/3. 

The d — 6 four fermion operator obtained from W^ 5 by higgsino dressing and involv- 
ing the neutrino has the form 

Oi= R & = j^Lrp (tfW) {d?y) , (4.10) 

where 



u u u d 



with u 6paa = -L (M^VckmV^) s (M° ag Vt KM P*) . (4.11) 

VyVd r \ /CUT 

As t R accounts for short distance renormalization factor of the corresponding RRRR d=5 
operator. Here we give values of those, which will be needed for further calculations: 

A s , R (u c , ^ 4.44 , A s , R (u c , t c ) ~ 4.0 . (4.12) 

4.2 Nucleon decay 

The operators responsible for p — > V p K + decay are 



^abc 



, ( [Cii2p( u adb)(s c v P ) + C v l2lp {u a s h )(d c v p ) + TZ p ll2p {u c a d c b )(s c u p ) + ^ 21p (M c a S c 6 )(4^p)] • 

(4.13) 

From these expressions we can calculate the partial widths for nucleon decay: 



T(p ->■ u p K + ) = 



(ftrf^ + a„7^ 12p ) (l + ^-(£> + 3F)) } 



B 

2 



(4.14) 



Here an, Ph are hadronic matrix elements and at /i — 2 GeV scale are [38] ~ |/3f/| ~ 
0.012 GeV 3 , while the values of other parameters are m p = 0.94 GeV, m K = 0.494 GeV, 
f n = 0.131 GeV, m B = 1.15 GeV, D = 0.8, F = 0.47. The factor R L ~ 1.25 is a long 
distance renormalization factor. 

Note that, different from the minimal SUSY S77(5) model, in Eqs. (4.4) and (4.5) the 
unitary matrix V appears. This matrix, by proper selection of its mixing angles, allows 
us to suppress proton decay so as to bring the partial lifetime within experimental limits. 
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Before demonstrating this with numerical results, in order to get a better feeling, we 
present an analytic study to leading order in certain small parameters. To leading order, 
let us ignore (i.e., set to zero) the 2 — 3 and the 1 — 3 mixing angles in the CKM matrix 
and in the V matrix. Let us also take the limit m u ,md,m e — > 0. In this limit, we get 

Cv nv nv rtv r\ 

1211 — U 1213 ~~ U 1121 — U 1123 ~~ U • 

Similar results hold for the corresponding 1Z V amplitudes. Therefore 

T(p ->■ V e K + ) = T(p V T K + ) = . (4.16) 

Only T(p — > V^K + ) will be non-zero due to the non-zero elements C\ 2X2 and C\ X22 which 
are given by 

r i212 = C U22 ~ g\ (l(u,d) + 7(«,e)) A a s e^\ s \^m9 c (sm9 c e l ^ +5 ^ + \> 21 e^) . (4.17) 

Note that in the limit V 21 — > the expressions of Eq. (4.17) will coincide with those of 
minimal SUSY SU(5). Now, we can select the matrix element V 2 i in such a way that 
these coefficients vanish (or are suppressed): sin# c e^ 2+52 ) + V 21 e t(i>1 = 0, or 

\V 21 \=sm6 c , Arg(V 21 ) = tt + <\> 2 + 5 2 - <pi . (4.18) 

With this conditions satisfied we get T(p — > V^K + ) ~ and the decay p — > VK + will 
be eliminated. Note that the conditions in Eq. (4.18) are easily satisfied. This is true 
for the second relation because all phases entering there are free. As far as the condition 
\V 2 i\ = sin# c is concerned, from (A. 12), with t\s 2 < 5tfs 2 we have \V 21 \ ^ftfs 2 . With 
the selection tfs 2 ~ 4 we get | V21 j ~ 0.2 ps sin# c . 

With the inclusion of 1 — 3 and 2 — 3 mixings, and vn u ^ e 7^ 0, the expressions get more 
lengthy, making analytical treatment harder. Thus, in the following we proceed with a 
numerical study, demonstrating the possibility of proton lifetime suppression. 



(4.15) 



4.3 Exact numerical results 

Following Eq. (2.22) we choose 



— cos Of , 9 2 = arccos — - cos 9 l 2 ) , 6t = arccos — cos 0i\ . 

m d J V m M / ' \ m r J 

(4.19) 



3 The elements 7^ ! j / 212 , 7\Li 12 2 are suppressed strongly and can be ignored. 
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Table 1: Particle masses (in GeV) obtained by the input given in Eq. (4.20) in MSSM. 

Then there are only three independent angles. We treat Of, 0\ and 0\ as free parameters 
and select them in such a way as to suppress d — 5 proton decay rates adequately. We 
also have the free phases 5\^, 001,2, 4>i,2, which we vary so as to suppress proton decay rate. 

For soft SUSY breaking parameters we adopt supergravity-inspired spectrum. How- 
ever, we deviate from mSUGRA and allow for non-universality in the Higgs boson mass. 
This is implemented by taking the pseudoscalar Higgs mass Ma and p, as independent pa- 
rameters. At the GUT scale we take as input, inspired by the "natural SUSY" spectrum 
ofRef. [39], 

M = 3 TeV, Mi /2 = 568.3 GeV, A = -5 TeV, 

tan/3 = 7, /i = 150 GeV, M A = 1 TeV, (4.20) 

where M (M1/2) is the usual universal soft mass for chiral matter superfields (gauginos) 
at the GUT scale, A the common trilinear term, while the Higgs sector is not universal 
{Mjj d 7^ M5 2 ). The value of tan f3 given is at the weak scale, corresponding to tan /3 = 6.75 
at the GUT scale. The parameters are chosen so that the SUSY spectrum is lighter than 
approximately 3 TeV, which can be discovered at LHC. For numerical calculations we 
used the code SuSpect [40], through which we make sure that the lightest (SM like) Higgs 
mass is ~ 125 GeV. The spectrum (at weak scale) we get for the input of Eq. (4.20) is 
given in Table 1. These values will be used in the calculation of proton lifetime. 

One choice of the three free angles and phases giving adequate suppression of proton 
decay rate is: 

Of = 1.3433, 0\ = 1.016, l 3 = 0.10275, 
<p l = 5 1 = 0, 2 = 3.3065, 5 2 = 1.883, 

wi = 2.515 , co 2 = 1-748. (4.21) 

With these input values we obtain for the decay rate p — > VK + 
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/0.012GeV 3 \ 2 /1.25\ V M T V . 

In Table 2 we summarize the partial lifetimes for this and other decay modes. Not all 
decay modes (induced by the d = 5 operators) are listed, those with lifetimes exceeding 
~ 5 • 10 36 years are not shown. Note that with further tuning of parameters, we may 
suppress even more the p — > uK + decay. However, we can not decrease much further the 
value of Mr because that would decrease the lifetime T~^ 5 (p — > fi + K°) whose value is 
already near at the experimental limit [41] (see Table 2). 

Note that with the value M T = 4.8 • 10 16 GeV (used in Eq. (4.22)), the mass of the 
£77(5) gauge bosons (X, Y) should be greater than about 2 x 10 16 GeV in order to be 
consistent with perturbtativity [26]. Such a value for Mx would mean that there is some 
chance for the observation of the gauge boson mediated nucleon decay such as p — > e + 7r°, 
but this will be challenging. 

One can try to increase the color triplet mass to further suppress the rates for the 
d = 5 modes. Due to the perturbativity constraint (see Eq. (2.11)) one needs first to 
increase the heavy gauge boson mass. For = m 8 this equals 

M x = M x /r 1/3 (4.23) 

where M x 2.10 16 GeV. By choosing r ps 1/10 or so M x and thus M T can be increased 
by a factor of 2. The color triplet mass can now be raised to M T ps 10 17 GeV, which would 
imply the scaling of all lifetimes for all modes in Table 2 upward by a factor of 4. Further 
increase of the triplet mass could jeopardize the expansion in inverse powers of the Planck 
scale, so we will not consider it. We see that, with the assumption that SUSY particles 
masses lie below about 3 TeV, which is testable at the LHC, proton lifetime cannot exceed 
about 2 x 10 34 years. This is within reach of ongoing and proposed experiments. 

We have not included gluino dressing of the effective d = 5 operators in order to 
obtain four fermion operators for proton decay. When universality is assumed, as we do, 
for the masses of the superpartners of the chiral fermions, the gluino dressing diagrams 
are highly suppressed [42] compared to the Wino dressing diagrams. This is primarily due 
to the antisymmetric nature of the QQQL operator in flavor. With the SUSY particle 
masses taken to be less than about 3 TeV, universality in the soft scalar masses is almost 
a necessity in order to suppress flavor changing neutral currents (FCNC) arising from the 
exchange of SUSY particles. If the third family squark and slepton masses are taken to 
be different from those of the (degenerate) first two families, FCNC processes may not be 
excessive. In this case, the gluino dressing contributions to nucleon decay may become 
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Table 2: Inverse widths for nucleoli decay. Calculations are carried out for the SUSY 
parameters (spectrum) given in Eq. (4.20), Table I. The model parameters are given 
in Eqs. (4.21), (4.19), along with M T = 4.8 x 10 16 GeV. Other parameters used can be 
found right after Eq. (4.14). 



r^ 5 (p 


4 • 10 33 yrs. 


r^(n uK°) 


2 • 10 33 yrs. 




1.0 • 10 34 yrs. 


r^b^^vr ) 


1.8 • 10 34 yrs. 




7.3 • 10 33 yrs. 


r d l 5 (n ^°) 


1.5 • 10 34 yrs. 



important, but typically the amplitude is not much more than that arising from the Wino 
dressing, see for eg. discussions in Ref. [43]. Thus, variation of SUSY spectrum would 
not significantly alter the upper limit on nucleon lifetime derived above, as long as the 
sparticle masses lie below 3 TeV or so. 

5 Conclusions 

In this paper we have shown that the main problems of the minimal renormalizable model 
based on SUSY SU(5) can be cured by adding a vector-like pair of 5 + 5 matter fields. 
This allows for the mixing of chiral families with the vector-like fields, which we show 
corrects the wrong mass relations of minimal SU(5). The mass splitting between the 
color triplets and the weak doublets of this vector-like fields improves the unification of 
the three gauge couplings. The color triplets from the 5# + 5# fields, which mediate 
d = 5 proton decay can have GUT scale masses, thus avoiding the rapid proton decay 
problem of the minimal model. The small number of couplings of this model enables us 
to make quantitative predictions for partial lifetimes for proton decay. We find that, in 
the favorable case that the LHC is sensitive to the discovery of the whole SUSY spectrum 
(corresponding to all the super-partner masses and Higgs boson masses < 3 TeV), at least 
some of the modes should have partial lifetimes shorter than about 2 x 10 34 yrs, which is 
within reach of proposed experiments. 
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A Deriving W e ff 

In this Appendix, we give details of obtaining the effective superpotential, both for the 
light fermion mass matrices, and for the d — 5 baryon number violating superpotential 
couplings. The effective superpotential is obtained by decoupling the extra heavy vector 
like states. First we integrate out the extra matter states. This is performed by block- 
diagonalization of the first two coupling matrices in Eq. (2.14). 

A.l Derivation of W mass 

With the transformation 



L = P X V}U , E c = P e cE c ' , D c = P dc V}D c ' , D = P q D' , 



(A.l) 



the matrices M? xA and M 4x4 get transformed to [30] 




(A.2) 




(A.3) 



The matrices in Eq. (A.l) are given by 




P 



e _<w 'Diag (e 




e 
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= e^ c Diag (e"^i d , <f*"* , e"^ , l) 
P q = Diag ( e i( *«rW , e * ( ^~ W , , l) (A.4) 



Z,rf 



c x —s{ 

e,d e,d e,d r> e,d e,d 

— S 1 S 2 c 2 U — C l S 2 

e,d e,d e,d e,d e,d e,d e,d e,d e,d 

c 2 s l s 3 ~ S 2 s 3 C 3 ~ C 1 C 2 s 3 



(A.5) 



i e,d e,d e,d e,d e,d e >d „ e >d e,d e,d t 

\ c 2 c 3 S l C 3 S 2 s 3 C l C 2 c 3 / 

where definitions for the entries of Eq. (A.5) see Eq. (2.20). We use the notation <px 
to denote the phase of a complex parameter X. Thus (f> yiVd is the argument of yiVj, etc. 
With all these, one can easily check that the matrices Me, and masses M^,Mc are 
given by Eqs. (2.19) and (3.1) respectively. The entries 0(vd) in Eqs. (A. 2), (A. 3) can 
be safely ignored. Thus, the diagonal block-entries in these matrices, together with Mjj, 
coincide with the terms of Eq. (2.18). 

A. 2 Deriving effective d = 5 operators 

Now we turn to the derivation of the effective d — 5 baryon number violating superpo- 
tential couplings. With the transformations of Eq. (A.l) and with 

q = P' q q' , u c = Py , (A.6) 

where 

P ^ = e -^ c Diag ^ M r^" d) , e^ M 2~^ 2 " d) , e^ M i ~^ 3 " d) ^ , (A.7) 
one can derive the couplings of the light states with the color triplets T, T: 

—l T M E P'qT + —u T MfjdT + — d cT M T D u c T + —e cT P'*M uU c T , (A.8) 
Vd v v ''./ v u 

where we have omitted primes for the quark and lepton states. The matrix P', without 

loss of generality, can be parameterized as: 

P' = Diag(e i5l ,e i52 ,l) . (A. 9) 

Further, integrating out the states T,T with mass M T) from Eq. (A.8) we derive the 
effective d = 5 operators given in Eqs. (4.2), (4.3). These are written in a flavor basis. 
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Finally, we present the matrices which appear in the d = 5 couplings written in the 

the mass eigenstate basis, using the transformations given in Eq. (2.24). These are the 
phase matrix P 

P = Diag(e itJl ,e iaJ2 ,l) , (A.IO) 

and the matrix 

V = VP , with V = V E P'Ul , P = Diag (e^ 1 , e i<t>2 , l) . (A. 11) 

The elements of the matrix V are: 

Vh ~ , V 12 ~ -^s'e* 1 + ^t\s e 2 e^ , 

m s rrifj. 



V ~ md +dddi5 1 m s m e eed diS 2 , m e,e S 3 

Vl3 — h c 2 s 3 e t l$2 t 2 s 3 e ~l > 

777,6 mb m r C\ 



V 21 ~ ^^e* - ^ a «e* , V 22 ~ , V 23 ~ -—t d 2 s d 3 e^ + ^t e 2 s e 3 , 
m s rrifj, mt, m T 

V„ ~ H^fl _ ^^.e e f d d iSi _ ^.e e e ii! 
'31 — n j t 2 6 3 t l i '2 e t i*3 t '2 e: ) 

m& c 2 m s m T m T 
V 32 ~ - ^s§e* , V 33 ^ 1 . (A.12) 

A. 3 An alternative derivation of W mass 

Here we provide an alternative, perhaps more intuitive, derivation of the effective mass 
matrices for the down-type quarks and charged leptons that follow from Eq. (2.16). We 
write down these matrices in a unified SU (5) notation, 

£ =(5. *.)(< ■ gift) (A.13) 



where 



m% = vMIh) (A. 14) 

M a = fi a + r] a (24 H ) , a = 1...4 (A.15) 

Here (24#) = 2v for the color triplet quark fields, while (2A H ) = —3v for the SU(2) L 
doublet lepton fields from the 5 a + 5 4 . Now we make a unitary rotation parametrized by 

(5i 5 4 ) ->• (5i 5 4 ) C/ (A.16) 
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U = U A 1 < A - 17 ' 



with 

A -Ax 

x T = (M U M 2: M 3 )/M 4 (A. 18) 

A = il + xx^)- 1 ' 2 , A = (1 + ^ x )-i/2 = (l + l^l 2 )- 1 / 2 (A.19) 

Note that the unitary matrix U is different for the quarks and leptons, since the Mj factors 
that enter into U are different. Similarly, the x, factors are not the same in these two 
sectors. We shall not explicitly show here the dependence of U or Xj on the fermion flavor, 
but it is to be understood. 

With the rotation of Eq. (A. 16) , Eq. (A. 13) becomes 

^{(x^Am% x^AM + AmJ (s^) (A ' 20) 

The heavy pair is now 5 4 — 5 4 , and the light mass matrices for down quarks and charged 
leptons become 

M D = A d m° M E = m°A eT (A.21) 

with 

Ui + 2riiV Lii — 3riiV , . 

x Di = ^ '±— , xeh = — — (A.22) 

/i4 + 2r)4V fJ,4 — 3l]4V 

where we have explicitly shown the separate matrices for down type quarks and charged 
leptons, using the GUT scale VEV v given in Eq. (2.2). 

The matrix A from (A.19) (for each sector separately) can be written explicitly as 

xx* 

A = 1- 



y/l+ |x| 2 (Vl + |x| 2 + l) 
'\ — CX\X\ —CX\X* 2 —CX\X\ 

-cx 2 xl I—CX2X2 — 0x2X3 I (A. 23) 

-CX3XJ — CX3X2 1 — CX3X3, 



with 



c = , v (A.24) 

y/l + |X| 2 [y/l + |X| 2 + lj 

The down quark and charged lepton mass matrices of Eq. (A.21) can be diagonalized 
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readily. Their eigenvalues are given by: 



2 i 2 i 2 

m 1 + m 2 + m 3 



2 2 i 2 2 i 2 2 

m 1 m 2 + m l m 3 + m 2 m 3 



N 2 (i + N 2 + N 2 ) + N 2 (i + N 2 + N 2 ) + N 2 (i + N 2 + |x 2 | 2 ) 

1 + \x\ 2 

W\d 2 \ 2 {\ + N 2 ) + |rf 2 | 2 M 3 | 2 (i + \xi\ 2 ) + |4| 2 Mi| 2 (i + N 2 ) 



1 + \x\ 



2 2 2 | 1 | ! 2 1 1 2 / a ot -\ 

m x m 2 m 3 = — - - — , (A.25) 

where dj's are common for Md and Me, while the different. From Eq. (A.25), 

it follows that realizing the mass hierarchy is possible only when \di\ are hierarchical, 
\di\ <^ | ^2 1 ^ |^3 1, m which case we can write down very simple formulas for the three 
masses: 

rrii = \di\ cos6>j . (A. 26) 

Here we define three mixing angles as: 

tan#! = tan 9 2 = — , ^ ^ :, tan# 3 = ! 3 ^ = (A. 27) 

Vl + ki| 2 V 1 + l^il 2 + l^al 2 

with < 9i < 7r/2. These are the same definitions used in Eq. (2.20). 

Noting that the mass matrix elements of Eq. (A. 13) can be all made real by redefini- 
tions of fields, we also obtain the unitary matrices that diagonalize M D and M E : 

U T M D V = M£ ag (A.28) 
V T M E U = Mg ag (A.29) 

We interchanged the notation [/■<->■ V passing from D to E, because it is that has 
the same form as M D . Again, the matrices U, V are different for down type quarks and 
charged leptons, we use the same symbol however. The unitary matrices U and V are 
given as (with \di\ <C \d 2 \ <C \d 3 \ 

( 1 -^1*2 -^t lC2 S 3 \ 

U - 1 ~^t 2 s 3 , (A.30) 

\mihS3. rn2 t2 s 3 1 / 

\ ?Tt3 C2 TU3 * ° / 

Y /ci + C 2 C 3 S1S2 -SlC 2 S 3 \ 

V ~ S1S2C3 c 2 + C3C1 -s 2 s 3 . (A. 31) 

1 + C1C2C3 \ . / 

\ S1S3 C1S2S3 C3 + C1C2/ 

Here q = cos^, Sj = sin6>j, tj = tan6>j. Terms of order (m^/m 2 ) and {m\/m 2 ) are ignored 
in the derivation of these matrices. 

It is possible to fit all quark and lepton masses consistently to the observed values. 
The mixing angles are related by the ratios given in Eq. (2.22). 
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